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ABSTRACT 

Context. Gamma-ray binaries have been established as a new class of sources of very high energy (VHE, >100 GeV) photons. These 
binaries are composed of a massive star and a compact object. The gamma-rays are probably produced by inverse Compton scattering 
of the stellar light by VHE electrons accelerated in the vicinity of the compact object. The VHE emission from LS 5039 displays an 
orbital modulation. 

Aims. The inverse Compton spectrum depends on the angle between the incoming and outgoing photon in the rest frame of the elec- 
tron. Since the angle at which an observer sees the star and electrons changes with the orbit, a phase dependence of the spectrum is 
expected. 

Methods. A procedure to compute anisotropic inverse Compton emission is explained and applied to the case of LS 5039. The spec- 
trum is calculated assuming the continuous injection of electrons close to the compact object: the shape of the steady-state distribution 
depends on the injected power-law and on the magnetic field intensity. 

Results. Compared to the isotropic approximation, anisotropic scattering produces harder and fainter emission at inferior conjunc- 
tion, crucially at a time when attenuation due to pair production of the VHE gamma-rays on star light is minimum. The computed 
lightcurve and spectra are very good fits to the HESS and EGRET observations, except at phases of maximum attenuation where pair 
cascade emission may be significant for HESS. Detailed predictions are made for a modulation in the GLAST energy range. The 
magnetic field intensity at periastron is 0.8±0.2 G. 

Conclusions. The anisotropy in inverse Compton scattering plays a major role in LS 5039. A simple model reproduces the observa- 
tions, constraining the magnetic field intensity and injection spectrum. The comparison with observations, the derived magnetic field 
intensity, injection energy and slope suggest emission from a rotation-powered pulsar wind nebula. These results confirm gamma-ray 
binaries as promising sources to study the environment of pulsars on small scales. 

Key words, radiation mechanisms: non-thermal — stars: individual (LS 5039) — gamma rays: theory — X-rays: binaries 



1. Introduction 

Gamma-ray binaries have been established in the past cou- 
ple of years as a new class of sources of very high energy 
(VHE, >100 GeV) photons. They are characterized by a large 
gamma-ray luminosity above an MeV, at the level of or ex- 
ceeding their X-ray luminosity. At present, all three such sys- 
tems known (LS 5039, PSR B1259-63 and LSI +61°303, re- 
cently possibly joined b y Cyg X-l) comprise a m assive star 
(lAharonian et alJl2005albl: lAlbertet alJl200ol 120071) . The com- 
pact object in PSR B1259-63 is a 48-ms, young radio pulsar. 
The VHE emission arises from the interaction of the relativistic 
wind from this pulsar, extracting rotational energy f rom the neu- 
tron s tar, with the stellar wind from its companion dTavani et alJ 
119941) . Particles gain energy at the shock bet ween the winds, re- 
sultin g in a small-scale pulsar wind nebula ( Maraschi & Treves! 
Il98ll) . The particles radiate away their energy as they are en- 
trained in the shocked flow, forming a comet-like trail of emis- 
sion behind the pulsar (Dubus 2006b|). 

The nature of the compact object and origin of the VHE 
emission remains controversial in LS 5039 and LSI +61°303, 
although recent observations indicate the radio emission of 
LSI +61° 303 behaves like the comet tail expected in the pulsar 
scenario dDhawan et al.l l2006). Alternatively, the VHE emission 
could originate from particles accelerated in a relativistic jet, the 
energy source being accretion onto a black hole or neutron star 



dDermer & Bottchetj l2006t iParedes et alJ 120061) . The rationale 
being that there is evidence for particle acceleration in the jets of 
microquasars and active galactic nuclei. However, hard evidence 
for accretion occuring in either LS 5039 or LSI +61°303 has 
been hard to come by (e.g. Martoc chia et al.ll2.Q05l) and the simi- 
larities between the three systems (and differences with the usual 
microquas ars) do not a rgue in favour of the accretion/ejection 
scenario (Dubus 2006q). 

Regardless of the actual powering mechanism, some par- 
ticles must be accelerated to high energies to generate the 
VHE gamma-rays. If these particles are leptons, the only viable 
gamma-ray radiation mechanism is inverse Compton scattering 
on the stellar photons. The massive stars in gamma-ray bina- 
ries have effective temperatures of several tens of thousand K 
and radii of about 10 R G , yielding luminosities of the order of 
10 39 erg s _1 . This provides a huge density of stellar photons in 
the UV band that VHE leptons may up-scatter, much greater than 
any other possible source of target photons (e.g. synchrotron or 
bremsstrahlung emission). 

The emitted VHE photons also have enough energy to pro- 
duce e + e~ pairs with the UV stellar photons. Most of the VHE 
flux may therefore be lost to the observer if the source is behind 
the star and VHE photons have to travel through the stellar light. 
Gamma-ray attenuation has been shown to lead to a modulation 
of the VHE flux with minimum absorption (maximum) at in- 
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ferior (superior) conjunction (Bo ttcher & Dermerl 120051: |t>ubus 
l2006ah . 

HESS observations have indeed shown a stable modulation 
of the VHE flux from LS 5039 on the orbital period with a 
maximum around inferior conjunction. This suggests attenua- 
tion plays a role and that the source of VHE gamma-rays cannot 
be more than about an AU from the binary (or attenuation would 
be too weak to modulate the flux). However, a non-zero flux is 
detected at superior conjunction where a large attenuation is ex- 
pected, possibly because of pair cascading. Moreover, the spec- 
tral changes that are reported do not fit with an interpretation 
based on pure at tenuation of a constant VHE source spectrum 
dAharonian et al.ll2006l) . 

Inverse Compton scattering also has a well-known depen- 
dence on the angle between incoming and outgoing photon. 
The photon flux from the star being anisotropic, the resulting in- 
verse Compton emission will depend on the angle at which it 
is observed. Hence, a phase-dependent VHE spectrum will be 
observed even if the distribution of particles is isotropic and re- 
mains constant throughout. Thi s effect has previou sly been in- 
vestigated in PSR B 1259-63 bv lBall & Kirk] | |2000l) who calcu- 
lated the radiative drag on the unshocked pulsar wind from scat - 
tering of stellar light, using results from Ho & Epstein (1989). 
The drag produces a Compton gamma-ray line with a strong de- 
pendence on viewing angle. 

This work purports to explain the HESS observations of 
LS 5039 using a combination of anisotropic inverse Compton 
scattering and attenuation in the simplest way possible. The 
aim is to constrain the underlying particle distribution and/or 
powering mechanism. §2 derives the main equations governing 
anisotropic Compton scattering in the context of gamma-ray bi- 
naries and discusses the principal characteristics to expect. §3 
presents the application to the case of LS 5039. The lightcurve 
and spectra observed by the HESS collaboration are reproduced 
by a model taking into account the photon field anisotropy and 
the attenuation due to pair creation. §4 concludes on the origin 
of the VHE emission from this system. 



2. Anisotropic Compton scattering 

Quantities in the electron rest frame are primed and quantities 
in the observer frame are left unprimed. The electron energy is 
y e m e c 2 , the energy of the incoming (stellar) photon is eo and the 
outgoing photon energy is e\ . These quantities are related in the 
electron rest frame by the standard 



1 + -2,(1 - COS©') 



(1) 



with 0' the angle between the incoming and outgoing photons. 
The incoming and outgoing photon energies are equal e[ = 
in the Thomson scattering approximation e' Q <k m e c 2 , or eo <c 
#*eC 2 /[y e (l - /3cosGq)] when expressed in the observer frame 
(00 is the photon angle with respect to the electron direction of 
motion). Scattering is also Thomson-like even if y e eo > m e c 2 
when the incoming and outgoing photon have almost the same 
direction (0' <K (2m e c 2 / 'e'^) 1 ^ 2 ). In the observer frame there is 
also an angle G CI h below which scattering will be Thomson-like. 
This angle is defined by 



1 / nipC" 
cos crit £ - 1 



y e £o 



(2) 



i.e. 6> cr it 5= 60° for yeo = 1 Me V. The cross-section in the electron 
rest frame is 

where r e is the classical electron radius and the photon energies 
e' Q j are related through Eq. ([TJ. 



2. 1 . Monoenergetic beam 

It is worthwile to consider first the simple case of a monoener- 
getic beam of photons scattering off a single electron. The main 
steps are listed below and a detailed derivation may be found in 
lFargionetal](ll997l) . 

In the observer frame, the incoming photon density (in 
sr^'erg" 1 ), normalised to the (constant) total photon density no 
(in photons cm -3 ), is 



dn 
dedQ, 



S(e- eo) 6 (cos G - cos Gq) 6(<p - <po) 



(4) 



with 5 the Dirac function. The frame origin is at the location 
of the electron (the frame orientation is arbitrary). The photon 
density in the electron frame is fou nd by using the invariance of 
dnlde dBlumenthal & Gouidll 19701) . 

The fraction of photons scattered per unit time , energy and 
solid angle in the electron frame is then given by dJones||1968; 
Blument haT& Gouidll 19701) 



dN' rC do- dn' , , 
= I I c dil de 

dt'de[dn\ JJ de[dQ.\ de'dQ! 



(5) 



which can be transformed to the observer frame using that the 
number of photons is invariant 



dN 



dN' dt' de[ d£l\ 
dtde\d£l\ dt'de'^dQ.^ dt de\ dCl\ 



(6) 



Q'j denotes the solid angle into which the outgoing photon is 
emitted and cos©' = cos G' cos G^ + sin G' sin ff x cos(</>j - <p'). 
Defining the polar angles 0o,i with respect to the direction of 
electron motion, the resulting differential photon spectrum is a 
function of y e , 6>o, 4>q, eo, 6\, <p\ and e\. The integration gives a 
rather unwieldy expression that can be found in the Appendix 
(Eq.Q. 

In the Thomson regime (e^ <s m e c 2 ), the outgoing photon 
energy is unequivocally related to the incoming photon energy 
since ej = e' Q . To each polar angle G\ corresponds a unique pho- 
ton energy. In the general regime there is also a dependence on 
the azimuth (see Appendix). Staying in the Thomson regime, the 
total spectrum emitted by an el ectron follows fro m the integra- 
tion over d£l' of Eq. (|5]l and is (Fargio n et al.ll 1997b 



nr 2 c 



dN 



dtde x Ipylen 



Te<, 



(7) 



where fii = cos6> and e\ varies between eo(l -y6^o)/(l ± fl)- 
This expression shows how the emitted spectrum depends upon 
the angle Gq between the monochromatic point source and the 
direction of motion of the electron. A more general expression 
is given in the Appendix (Eq. IA.6l l. 
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Fig. 1. Dependence of the viewing angle on the inverse Compton spectrum. The source of photons is a star with kT-l eV. The 
electron cloud is situated at a distance d - 2/?*. The electrons are distributed according to a power-law dn e - y~ 2 dy e . The left 
panel shows the variation of the spectrum with angle when the interaction occurs in the Thomson regime (electron energy range 
10 3 < y e < 10 5 ). In the right panel the interaction occurs in the Klein-Nishina regime (electron range 10 5 < y e < 10 7 ). In each 
panel, the spectrum is shown at viewing angles if/ =15° (bottom), 30°, 60°, 90°, 120° and 180° (top). The observer sees the electron 
cloud in front of the seed photon source when the angle is small. Solid lines are calculated taking into account the finite star size 
(Eq.[T3l; dashed lines correspond to the point source approximation (Eq.[T2l. 



2.2. Kernel for spectral calculations 

The monochromatic, single photon result can be used as a kernel 
to integrate over general electron and incoming photon distribu- 
tions. The total spectrum in photons s _1 erg -1 sr~' is then given 
by 



dN« 



dtdei c/Q 



'A 



Iff! 



dN 



dtdei dQ. 



nodClodeo 



dy e dQ., 



■dQ. e dy e (8) 



where the evaluation of the kernel must take into account the 
changes in electron direction with respect to the given direction. 
However, this expression can be simplified. 

The electron energy must be very large y e » 1 in order to 
emit VHE photons. The emission is strongly forward boosted in 
the direction of the electron motion by relativistic aberrations. 
An observer looking at the inverse Compton emission from an 
isotropic cloud of relativistic electrons sees essentially only the 
emission emitted by those electr ons moving within an angle l/y e 
from the line-of-sight (see e.g. iBall & Kirk 2000). Their emis- 
sion is almost entirely focused into the line-of-sight. Photons 
emitted slightly away from the line-of-sight and included in the 
integration compensate to order l/y e for the emission from elec- 
trons moving at larger angles. Therefore, to a good approxima- 
tion, 



dN 



dn. 



-d£l e 



dN 



dtde\dD.\ dy e d£l e dtde\ 
and the spectrum will be given by 



drip 



dN tc 



dtde\d£l\ 
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dN 



dtdei 



a dy e dQ. t 



dn e 



(9) 



,dy e dQ e 



dy e 



where the kernel is given by Eq. (0 or Eq. ( |A.6t , evaluated at 
the angle a between the point-like photon source, the electron 
cloud and the observer. If e b s is a unit vector from the electron 
cloud to the observer and eo is a unit vector from the photon 



source to the electron cloud, expressed using 6q and <po, then 
fi a = cos a = eo.e obs . 

For scattering on an isotropic distribution of photons, e bs 
can be arbitrarily oriented so that a = 6q. For a blackbody of 
temperature T*, 



nodQ.Q = 



/z 3 c 3 exp(«b/*r*)-l 



(11) 



For a point-like star of radius at a distance d+ from the elec- 
trons, e bs can be defined on the plane containing the three loca- 
tions so that, again, a = 9q. The photon distribution is 



lh,(IQu = /I) ^ | ,/0()<//„ - //,:', )<)((/)„ ).-/o<i 



(12) 



with fo as defined in the previous equation and where if/ is the 
angle between the star centre, the cloud and the observer. The 
integral on Qq is direct so the kernel only needs be numerically 
integrated on £q and y e . Finally, for a star of finite size, the inte- 
gration element is 

n dQo = fo cos 6 dQo, (p e [0,2n], sm6 e [0, (13) 

and fi a = cosi^cos^o + sini^sinyuocos^o- This requires a 
quadruple numerical integral. 

The electron distribution will be assumed to be isotropic in 
the following so that the expression in Eq. (0 is a function f e 
of y e only and J f e dy e gives the total number of electrons per 
steradian. 



(10) 2.3. Anisotropic scattering of stellar photons 



Figure 1 shows example calculations of the inverse Compton 
spectrum from a distribution of electrons scattering photons 
emitted by a star, as seen from different viewing angles. The 
incoming photons have a blackbody distribution and the elec- 
trons have a power-law distribution dn e = y~ 2 dy e . The viewing 
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angle if/ is defined as the angle between the star, electron cloud 
and observer (cos if/ = e*.e b s with e* a unit vector from the star 
centre to the cloud). Two cases are shown: one corresponding to 
scattering in the Thomson regime and one for the Klein-Nishina 
regime. For each case, results obtained in the point source ap- 
proximation and taking into account the finite size of the star are 
compared. 

When scattering occurs in the Thomson regime (e ( , <K m e c 2 ), 
the maximum energy j 2 €q(\ + 1 -/Jcos if/) decreases with de- 
creasing viewing angle if/, i.e. when the electrons move in front 
of the star as seen by the observer (left panel of Fig. [TJ. This is 
to be expected as the electrons are then forward scattering radi- 
ation that is less energetic in their rest-frame than in the head-on 
case because of the 1-/3 cos iff term in the Lorentz transform. 
The other effect is a lower rate of emission for low if/ (as can be 
directly deduced from Eq.|6]and seen in the left panel of Fig. [T}. 
This is also due to the decrease in the density of incoming pho- 
tons in the electron rest frame when both particles move in the 
same direction. Scattering is more likely to occur wh en the par- 
ticles collide head-on (e.g. lSazonov & S unvaev 2000). 

These effects are pronounced in the point source approxima- 
tion and are diluted when taking into account the finite size of 
the star (see dashed lines compared to full lines in Fig.[TJ. With a 
star of finite size, electrons see incoming photons from a variety 
of angles, which contributes to raising the seed photon density 
in the electron rest frame when if/ = (and to slightly decreas- 
ing it at = n). Because the density is tied to 1 - /? cos if/, this 
suggests a simple rule-of-thumb, corroborated by numerical in- 
vestigations: the effect of the finite star-size should be taken into 
account when sin if/ ^ R^/d* but can otherwise be neglected. If 
the observer is within the cone defined by the star with the elec- 
trons at apex, then the density of photons seen by the electrons 
moving towards the observer will be significantly greater than in 
the point source case. Outside of this cone, the difference with a 
point source approximation is minor. In Fig. [1] the star angular 
size seen by the electrons is 30° (defining the cone opening an- 
gle) and the point source approximation is indeed acceptable for 
if/ > 30°. 

When scattering occurs in the Klein-Nishina regime (e^ > 
m e c 2 ), the maximum energy is almost constant at y mw jn e c 2 re- 
gardless of viewing angle. For large viewing angles, the spec- 
trum is soft due to the decrease in cross-section in this regime, 
just as in the isotropic case. At small viewing angles, the seed 
photon energy in the rest frame of the electron is lower than 
in the head-on case because of the angle dependence in the 
relativistic boost, as described above for the Thomson regime. 
Moreover, since the limit between Thomson and Klein-Nishina 
regimes is at eoTe(l -yScos if/) as m e c 2 , scattering can reach back 
to the Thomson regime for small enough viewing angles, regard- 
less of the electron energy (see Eq.[2]in §2). There are two con- 
sequences. First, the amplitude of the variations with viewing 
angle is smaller than in the Thomson regime, because at small 
if/ the decrease in photon density is compensated by the larger 
cross-section. Second, since there is no drop in cross-section at 
small if/, there can be a significant hardening of the spectrum 
compared to the spectrum at larger if/ (right panel of Fig. Q]). 
These spectral effects may play an important role in modelling 
the emission from gamma-ray binaries, for which scattering oc- 
curs mostly in the Klein-Nishina regime. This is investigated in 
the next section. 



3. Application to LS 5039 

The influence of anisotropic scattering on the emission from 
gamma-ray binaries can be sketched from the results of the 
previous section. If the high energy emission is due to inverse 
Compton scattering off electrons co-rotating with the binary, the 
viewing angle of the observer will vary with orbital phase, induc- 
ing changes in the observed spectrum — all other things being 
set equal (particle distribution and location, distance to the star 
etc). 

Anisotropic scattering will most influence the emission from 
systems with high inclinations, if the electrons are located in the 
orbital plane. At low inclinations the changes are expected to be 
minor as the scattering angle if/ stays close to n/2. On the other 
hand, for high inclinations the inverse Compton spectrum may 
change significantly between inferior and superior conjunctions. 
The emission will be intense and soft at the time of maximum 
attenuation by pair production, and low and hard at the time 
of minimum attenuation. Anisotropic inverse Compton emission 
combined with attenuation of VHE photons can therefore play 
an important part in (1) reducing the amplitude of the varia- 
tions expected from a simple attenuation model; (2) hardening 
the spectrum at high flux states compared to expectations from a 
calculation assuming an isotropic flux. 

LS 5039 presents an ideal testbed. The massive star has an 
06.5V spectrum (T* = 39,000 K, fl*=9.3 R Q , M*=23 M ) in 
a 3. 9 day eccentric orbi t (e = 0.39) with its compact compan- 
ion dCasares et al.l l2005). A diagram of the binary orbit oriented 
on the sky is shown in Fig. [2] The measured radial velocity of 
the O star constrains the inclination to about 60° for a neutron 
star companion and about 20° for a black hole. The compact 
star moves from one to three stellar radii from the surface of the 
massive star. 

The intensity and spectral variations have been well- 
established in LS 5039 by HESS observations, concluding that 
pure attenuation of a constant VHE spectrum could not ex- 
plain the observations to satisfaction (see §1). Given the above 
discussion, this section examines whether taking into account 
anisotropic scattering provides an improved agreement. 

3.1. The radiating electrons 

Two main assumptions are made to calculate the emission. First, 
the electrons are assumed to scatter radiation at the location 
of the compact object, in a small region compared to the or- 
bital separation. This is a very good approximation in the pulsar 
wind nebula scenario where the highest energy electrons emit 
the gamma-ray radiation close to the shock. (The cooled elec- 
trons then emit in radio well away from the system.) This may 
or may not be appropriate in the case of a relativistic jet, where 
emission can occur at various distances along the outflow. This 
is further discussed jj3.4| 

Second, the adopted distribution of particles is the steady- 
state distribution for constant injection of particles, taking into 
account synchrotron and inverse Compton losses. The magnetic 
field in the radiating zone is assumed to be homogeneous. The 
radiative losses occur on very short timescales compared to the 
orbital timescale so the steady-state approximation is justified 
except for low energy particles whose radiative timescale be- 
comes longer than their escape timescale from the radiating 
zone. This occurs at y e ~ 10 3 (see below). The injection spec- 
trum is a power-law dn e oc y e p dy e with an exponential cutoff at 
the maximum y max allowed by comparing acceleration and ra- 
diative timescales. 
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Fig. 2. The binary orbit of LS 5039 as seen from directly above. 
The 06.5V star radius is to scale. The binary orientation is set 
for an observer at the bottom of the diagram. The binary inclina- 
tion on the plane of the sky is not taken into account. The num- 
bers indicate the orbital phase (mean anomaly) at various posi- 
tions. Periastron passage is indicated by a full line (orbital phase 
0orb=O). The dashed line is the line of conjunction s ((p SUD » 0.06, 
0i n f a 0.72). The orbital parameters are taken from lCasares et all 
(2005). 
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Fig. 3. Steady-state electron distribution N s along the orbit of 
LS 5039. The injection spectrum is a yT 2 power-law with an 
exponential cutoff at y e » 10 8 (see §3. 11 1. The magnetic field 
varies as B — 0.8c/ j G, where c/o.i is the orbital separation 
in units of 0.1 AU. Inverse Compton losses in the Thomson 
regime (y e < 7kn) and synchrotron losses (y e > ys) steepen 
the index of the injected distribution by one to N e «- yT 3 ■ 
Inverse Compton losses in the Klein-Nishina regime dominate 
between 7kn < Te < y s, causing a hardening of the distribution 
dModerski et al.ll2005l) . The steady-state distribution varies little 
with orbital phase since ys 00 (fidy 1 stays constant: the changes 
with orbital phase produce only a slight thickening of the line in 
the above figure. 



The minimum acceleration timescale for TeV electrons 
(y6=10 6 ) is set by the gyrofrequency and is f acc * 0.06 ye/B\ s 
with B\- \ G the magnetic field intensity. The synchrotron cool- 
ing timescale is ?s ~ 770/Bjye s. For electrons with Lorentz 
factors y e > yKN ~ 6 10 4 r + ' 4 , inverse Compton scattering of 
stellar photons occurs in the Klein-Nishina regime. The corre- 
sponding timescale is he ~ 2 y^d^ x / [In yg + 1 .4] (TV^-^io) 2 s 
dBlumenthal & Gouldfl970l) with f M =40,000 K, tf*'io=10 R e 
and do. i is the orbital separation in units of 0. 1 AU (the LS 5039 
orbital separation at periastron). 

The steady-state distribution derives from a comparison of 
these three timescales. Synchrotron losses dominate over inverse 
Compton losses above a critical ys given by (fs=fic): 

ys « 6 ■ 10 6 (r*, 4 /?*,io)/CBi*.i). (14) 

At the highest energies, y max is therefore set by synchrotron 
losses (f a cc=fsX which gives y max « 1.2 lO 8 /?^ 1 ^ 2 . Assuming 
continuous injection of electrons with a yT p spectrum, the 
steady-state distribution is steepened by synchrotron losses be- 
tween ys and y max to a y e p 1 power-law. Inefficient Klein- 
Nishina losses dominate between yKN and ys, producting a hard 
spectrum mirroring the decrease in energy loss rate with increas- 
ing y e in the Klein-Nishina regime. Below y^N inverse Compton 
losses in the Thomson regime result in a y"/^ power-law as in 
the synchrotron case. 

Steady-state distributions obtained using a full numerical 
calculation follow very w ell the main character istics outlined 
above (Fig. [3] see also iModerski et all 120051) . The inverse 
Compton losses are treated in the isotropic approximation since 
the magnetic field will quickly randomize particle directions. 
The particles see, on average, the equivalent of an isotropic radi- 



ation field; but the inverse Compton spectrum received by an ob- 
server at a fixed location changes with viewing angle. In Fig. [3] 
the injection is a power law y e p with p — 2 and the distribution 
between yKN and ys is roughly proportional to y~ 1 - 3 . The slope 
of this distribution depends on the slope of the injected spec- 
trum. For power-law injections y e p with hard indices (p < 2) the 
slope between y^ and ys tends to yT 1 . For soft indices p > 2, 
the hardening gradually disappears, reaching yT 2 between yKN 
and ys for an injection with p = 3. As discussed below, the ob- 
servations of LS 5039 constrain p to about 2. 

This steady-state distribution is a very good approxim ation 
to the more detailed pulsar wind model of Dubus (2006b]) for 
electrons with y e ^ 10 3 : lower energy electrons escape from 
the vicinity of the pulsar without radiating much of their energy. 
More generally, this distribution should apply equally well to 
any leptonic model assuming a constant injection of non-thermal 
particles cooling in the vicinity of the compact object via syn- 
chrotron and inverse Compton radiation. 

3.2. Compact pulsar wind nebula: orbital lightcurve 

With the inverse Compton losses fixed by the geometry, the only 
remaining free parameters are the slope of the injected power- 
law, the total energy in radiating electrons and the value of the 
magnetic field. In the case of a compact pulsar wind nebula, the 
magnetic field is determined by the conditions at the pulsar wind 
termination shock. Its intensity sets ys, which in turn will fix the 
frequency above which a break will be seen in the VHE gamma- 
ray spectrum. In principle, B may vary with orbital phase as the 
eccentric orbit brings the pulsar at various radii in the stellar 
wind. However, the magnetic field intensity is inversely propor- 
tional to the shock distance from the pulsar, and the latter is 
roughly proportional to the orbital separation so that B oc l/d 
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Fig. 4. Predicted orbital lightcurves for LS 5039 in the case of 
a neutron star (i = 60°). Top panel: the integrated photon rate 
above 1 TeV (full line) due to anisotropic inverse Compton scat- 
tering and the transmission exp(-Tyy) for the pair production 
process, also integrated above 1 TeV. Inverse Compton scattering 
is minimum at inferior conjunction (0 or b ~ 0.72, see Fig. [2j». The 
absorption due to pair production is also minimum at this time. 
Bottom panel: the resulting orbital lightcurve (full line) com- 
pared to the HESS observations. Combining anisotropic inverse 
Compton emission and attenuation by pair production produces 
a peak at or b « 0.8 consistent with the observations. The agree- 
ment is good except at periastron where cascade emission (ig- 
nored here) may be important. The dashed line shows the photon 
rate in GLAST above 1 GeV (ph crrT 2 s _1 ). The model predicts 
a peak in the GLAST lightcurve close to periastron and a min- 
imum at inferior conjunction. The normalizations are arbitrary. 
The lightcurves are calculated using the electron distributions 
shown in Fig. [3] 



(see e.g. Dubus 2006b). In this case, the distribution of particles 
will not change along the orbit as ys 00 {Bd)~ l . 

Figure [4] shows the expected lightcurve at different orbital 
phases with B = 0.8 G at periastron and p = 2 (using the elec- 
tron distributio n shown in Fi g. |3j. The orbital elements were 
computed as in lDubusI (l2006al) . The unabsorbed intensity is high 
close to superior conjuction and small at inferior conjunction, as 
explained in §2.31 The angle to the observer varies between 30° 
and 149° whereas the angular size of the star at the compact ob- 
ject is 30° at most: the finite size of the star, taken into account in 
the calculation, has a minor effect on the results. The attenuation 



lightcurve, computed following lDubusI d2006al) is also shown. It 
peaks at inferior conjunction where attenuation is minimum. 

The lightcurve including both anisotropic emission and at- 
tenuation by pair production reproduces very well the observed 
lightcurve. Most notably, the combination of low attenuation, in- 
creasing photon density and a hard inverse Compton spectrum 
produces a small peak after inferior conjunction that appears to 
be present in the HESS observations. The peak is a key feature of 
this model. This lightcurve is very robust against changes in the 
value of the magnetic field used, or even in the type of particle 
distribution used. At higher inclinations, a weaker peak appears 
before inferior conjunction as the variations in viewing angle 
cause a larger drop in inverse Compton emission at or b = 0.72. 
However, this model still predicts little to no flux at and after 
periastron because of the very strong attenuation of the emission 
emitted around the pulsar. A possible explanation is that a pair 
cascade develops. 

The lightcurve above 1 GeV is also plotted in Fig. [4] 
Attenuation is negligible and the variations mostly follow 
the photon density modulo some modifications due to the 
anisotropy: for instance, the minimum is at inferior conjuction. 
GLAST should therefore see a modulation in the flux from 
LS 5039 with a peak close to periastron and a minimum at infe- 
rior conjunction, almost anti-correlated with the HESS modula- 
tion. 

A similar lightcurve has been obtained by iBednarekl d2007l) . 
using a complex Monte-Carlo code simulating the effects 
of anisotr opic scatt e ring a nd the development of cascades. 
However, Bednarekj (120071) wrongly interpreted the GLAST 
modulation as being due to stronger cascade emission close to 
periastron. As described above, the modulation is due to a com- 
bination of increased seed photon density and anisotropic effects 
and not to cascade emissiorQ. 

3.3. Compact pulsar wind nebula: phase-resolved spectra 

Figure shows the evolution of the attenuated and unattenu- 
ated spectra with orbital phase. These were used to produce the 
lightcurves shown in Fig. [4] The spectra display a complex in- 
terplay between the varying threshold for pair production, the 
high absorption it causes at superior conjunction when the in- 
verse Compton flux is high, and the weaker but harder inverse 
Compton emission at inferior conjunction. The variations in the 
GeV (GLAST) range have a very large amplitude, with a flat 
spectrum at the highest intensities and a hard spectrum at low 
intensities. This should easily be accessible to GLAST in the 
very near future (Dubois 2006). Note that synchrotron emission 
contributes significantly to the emission below a GeV and that 
this is not taken into account in this lightcurve. Its impact is to 
soften the spectrum and reduce the amplitude of the variations 
below a GeV (see §3.4 below and Fig.|6]l. 

The attenuated spectrum averaged over the full orbit is 
shown in Fig. [6] The hitherto puzzling drop between the EGRET 
and HESS spectra is very well reproduced by the model with- 
out invoking a cascade. The inverse Compton spectrum by itself 
underestimates the EGRET flux by factors of a few but, taking 



1 Bednarek (2007) also confused the phases of inferior and superior 
conjunctions (Fig.[2j. The compact object is on the near side of the or- 
bit (inferior conjunction) at phases 0.4-0.8 so that the broad maximum 
is not due to the stronger Compton scattering expected when the ob- 
ject is behind the star (see Fig. [5}. Similarly, the dip at phase 0.7 is not 
due to stronger absorption (expected at superior conjunction) : it actu- 
ally occurs at the phase of minimum absorption and is due to the lower 
Compton emissivity at inferior conjunction, as described above. 
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Fig. 5. Evolution of the model inverse Compton spectrum with orbital phase in LS 5039 (neutron star case). The intrinsic emission 
spectra are shown with full lines and the dashed lines show the spectra after attenuation by pair production on stellar photons during 
the propagation of the gamma-rays through the system. The underlying electron distributions are those shown in Fig. [3] Left panel: 
from top to bottom the spectra correspond to orbital phases orb =O.O3, 0.09, 0. 15, 0.24, 0.34, 0.44, 0.56 and 0.66 (see Fig.|2|. Right 
panel: the plotted orbital phases from bottom to top are 0.66, 0.76, 0.85, 0.91, 0.97 and 0.03. 



into account the synchrotron emission from the electrons using 
the adopted magnetic field intensity (B=0.8 G at periastron and 
varying as l/d), the calculated synchrotron emission produces 
a very good match to both the EGRET and HESS spectra. Note 
that the average HESS spectrum is not shown for reasons of clar- 
ity in Fig.|6]but is close to the 'high' state spectrum (see below), 
with a slightly higher luminosity. 

The two average spectra for the phase intervals of the HESS 
'high' (0.45 < orb < 0.9) and 'low' state ( orb < 0.45 or 
(ftorb > 0.9) spectra are also shown in Fig. l6l ( Ahar onian et alj 
2006). Reproducing the cutoff in the high-state HESS spectrum 
strongly constrains the magnetic field intensity to ss 0.8 G at pe- 
riastron. A higher magnetic field moves the cutoff to lower ener- 
gies and is inconsistent with the data. A lower B moves the cutoff 
to higher energies and hardens the spectrum too much. The high- 
state spectrum is rather sensitive to the value of B: the acceptable 
range is only B = 0.8 + 0.2^' G. Outside of this range the fit 
does not go through the error bars of the HESS data points. 

The synchrotron emission contributes significantly below 
1 Ge V, diluting the hardening of the spectrum around (p = 0.7 ex- 
pected from pure inverse Compton emission. Actually, a soften- 
ing is predicted below a few GeV. The GLAST lightcurve shown 
in Fig. [4] is not noticeably changed (on a linear scale) by tak- 
ing synchrotron emission into account. The hard electron distri- 
bution, naturally resulting from inefficient Klein-Nishina losses 
here, is instrumental in obtaining the flat spectrum in the HESS 
range. The range 7kn < Te < Ts of this hard distribution de- 
pends upon the value of the magnetic field, but its shape is inde- 
pendently set by the index p of the injected power-law y e p . With 
p 5 1.7 the predicted HESS spectrum is too hard and the emis- 
sion in the EGRET band is too low. With p <: 2.3, the predicted 
HESS spectrum is too soft and the EGRET emission is too large. 
Therefore, the slope of the injected power-law is constrained to 
p = 2 + 0.3. 

Besides the magnetic field intensity and slope of injected 
electrons, the other free parameter is the normalization of the 
electron distribution. The fit was obtained for a total energy in 
electrons from y e = 10 3 to +oo of 3T 37 erg. This energy corre- 



sponds to the injection of 10 erg s in particles, assuming an 
escape timescale from the radiative zone of 30 s (longer than the 
radiative timescale under consideration). In the pulsar wind neb- 
ula the shocked electrons have a bulk velocity « c/3 so that the 
escape timescale corresponds to a radiating zone of 3 10 11 cm, 
comparable to the shock si ze found for typical wind parameters 
in LS 5039 (Dubus 2006b). The estimated injection energy rate 
is consistent with a reasonable p ulsar spindown power, s uch as 
that measured in PSR B1259-63 (Manchest er et al .1119951) . 



The low-state spectrum is responsible for most of the orbit- 
averaged emission in the EGRET range, which is nicely fit by 
the model. However, the HESS low-state spectrum is not satis- 
factory. This spectrum corresponds to phases where the intrinsic 
inverse Compton emission is both soft, as the observed spec- 
trum, and intense. The intrinsic emission is actually strongest at 
the times of highest attenuation so that the two effects compen- 
sate somewhat. However, the cross-section for pair production 
drops above a few TeV. Therefore, the predicted phase-averaged 
low-state is not a pure power-law but still shows hints of an atten- 
uation line with a kink at high energies. Changes in the electron 
distribution may also help to reduce the discrepancy. A cutoff at 
a lower y e (i-e. a higher magnetic field) than that shown in Fig. [3] 
would yield a better agreement if it occurred at the appropriate 
orbital phases. However, at this stage it appears more reasonable 
to investigate first the impact of pair cascading on this spectrum, 
as this is required to explain the detection at periastron. 

The model contains only three parameters: the slope of the 
injected power-law, the particle distribution normalization and 
the magnetic field intensity at periastron (or any other arbitrary 
orbital phase). The shape of the particle distribution and the asso- 
ciated emission along the orbit are then unequivocally predicted. 
The parameters were adjusted so as to fit the high-state HESS 
spectrum. That this choice also fits very well the EGRET obser- 
vations gives strong support to this simple-minded model, even 
if the low-state HESS spectrum is not reproduced to satisfaction. 
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Fig. 6. Comparison with the EGRET and HESS observations 
of the LS 5039 model spectra for a neutron star. The EGRET 
bowtie is in dark grey a nd the HESS high-state and low-state 
bowti es are in light grey (Hartm ali et al.ll 19991: lAharonian et al.l 
2006). The corresponding HESS deconvolved spectral points are 
also shown (with a dot identifying the low-state points). Fluxes 
have bee n transformed to lu minosities assuming a distance of 
2.5 kpc dCasares et alj |2005). The full grey line is the average 
spectrum calculated using the results of Fig. [5] It reproduces 
well the drop in flux from EGRET to HESS (the average HESS 
spectrum is close to the high-state spectrum shown). The high- 
state spectrum (full dark line) is very well reproduced provided 
the magnetic field at periastron is lower than 0.8 G. The low- 
state spectrum (dashed dark line) is not reproduced well, possi- 
bly because cascade emission contributes significantly at these 
orbital phases where pair production is very important or be- 
cause the electron distribution varies along the orbit. Here, the 
synchrotron emission from the electrons is taken into account 
with B = 0.8c/ | G as derived from the VHE spectrum. Its con- 
tribution to the spectra is shown by the dotted lines (from bottom 
to top: low-state, high-state, and orbital average). 



3.4. Black hole jet? 

This subsection examines how the results are changed if the 
compact object is a black hole. The main effect is that con- 
sistency with radial velocity curves require the inclination to 
change to 20° (4.5 M Q black hole). The variation in viewing an- 
gle is then reduced to the interval 70°— 110°. The electrons are 
still assumed to be accelerated in the vicinity of the black hole 
and to reach a steady-state distribution such as the one described 
above. Here, the magnetic field has a fixed value as there is no 
a priori reason for it to change with the orbital separation. This 
gives a moderate change of a factor 2 in the break ys of the elec- 
tron distribution, because the orbital separation varies by a factor 
2, in contrast to the situation described in Fig. [3] 

The orbital lightcurve and the spectra obtained with B = 
0.8 G and p = 2 are shown in Figs. |7][8] In contrast with the 
neutron star case, there is only one broad peak in the predicted 
HESS lightcurve. This is because the reduced variation of the 
viewing angle does not lead to a large drop in scattered flux at in- 
ferior conjunction. The small peak predicted at high inclinations 
(neutron star) can therefore be used as a discriminant between 



the two cases. The averaged spectra are much harder than in the 
neutron star case. The amplitude of the variation at GeV energies 
is less than for a neutron star and the average flux overestimates 
the EGRET emission. The poor fit of the low-state spectrum re- 
mains. Both the lightcurve and spectra are arguably not as good 
fits as those obtained in the neutron star case, but not so much as 
to exclude that LS 5039 is seen at a low inclination (and hence 
contains a black hole). 

Emission from a relativistic jet may differ from the estimate 
above. Any Doppler boosting will change the observed spec- 
trum. However, the resolved radio emission, if interpreted as a 
compact jet, implies only a moderate velocity and little boosting 
dParedes et alJ l2000). Modest Doppler (de)boosting may also be 
expected from the pulsar wind emission as its post-shock speed 
is approximately c/3. More importantly, emission may occur all 
along the jet and not just be localized near the black hole. Far 
from the compact object, the viewing angle tends to become the 
inclination angle (ifr — > i) regardless of orbital phaseQ. Hence, 
emission at progressively higher altitudes in the jet is less and 
less influenced by anisotropic effects. The emission is also less 
attenuated by pair production, with r yr negligible at heights 
^ 1 AU. If most of the emission occurs far in the jet, and assum- 
ing the electron distribution stays constant, the flux modulation 
is only linked to the stellar photon density. The result is a con- 
stant spectral shape, peak flux at periastron and a trough at apas- 
tron. These are inconsistent with the observations. Therefore, a 
jet model for LS 5039 probably requires either (1) that most of 
the emission occurs close to the compact object in order to re- 
produce the orbital gamma-ray modulation via anisotropic scat- 
tering and attenuation or (2) that the emission occurs away in 
the jet and that some unspecified intrinsic mechanism changes 
the particle distribution and/or the radiation process. 

4. Conclusion 

The anisotropic behaviour of inverse Compton scattering has a 
major influence on the emission from gamma-ray binaries. In 
these sources, the massive star provides a large source of seed 
photons with energies around an electron-volt. If high energy 
electrons are accelerated in the vicinity of the compact object, 
then the angle between the star, compact object and observer 
changes with orbital phase. The variation in viewing angle leads 
to a strong modulation in both the intensity and spectral shape 
of the scattered radiation. 

Scattering stellar photons to the TeV range requires very en- 
ergetic electrons with Lorentz factors y e ~ 10 6 - 10 7 . The scat- 
tering therefore occurs in the Klein-Nishina regime. In this case, 
the anisotropy results, at inferior conjunction, in a harder and 
fainter spectrum than predicted using an isotropic approximation 
for the incoming photons. Crucially, inferior conjunction also 
corresponds to the phase at which the produced VHE gamma- 
rays are less attenuated by pair production on stellar photons. At 
other phases the emitted spectrum is close to the one obtained us- 
ing the isotropic photon field approximation and can be severely 
attenuated by pair production. The result is a complex interplay 
that reduces the amplitude of the variations expected from a pure 
attenuation model and a hardening at inferior conjunction. 

2 Note that two errors have slipped by in Dubus ( 2006a) when dealing 
with the case of a VHE source perpendicular to the orbital plane. In the 
last equation of A. 2 the angle for emission perpendicular to the plane is 
given as cos if/ = (do Id) cos if/g = (do Id) sin sin but should be cos i/r = 
(do /d) sin 8 sin i - (z/d)cosi. The other is that Fig. 8 (attenuation with 
height) was calculated at a fixed viewing angle of 76° . The conclusions 
are unchanged. 
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Fig. 7. Predicted orbital lightcurves for LS 5039 in the case of 
a black hole (i = 20°). The full line is the integrated photon 
flux above 1 TeV (HESS), the dashed line is integrated above 
1 GeV (GLAST). The variations in viewing angle are reduced 
compared to the high inclination (neutron star) case (Fig.[4]i and 
there is only one broad maximum in the HESS lightcurve. The 
electron distribution is calculated as described in Fig.[3]but using 
a constant magnetic field intensity of 0.8 G. 
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Fig. 8. Comparison with the EGRET and HESS observations of 
the LS 5039 model spectra for a black hole. The radiating elec- 
trons are injected in the immediate vicinity of the black hole. 
The magnetic field intensity used to fit the high-state spectrum 
is 0.8 G, constant throughout the orbit. The injected electrons 
have a power-law of index p — 2. The line coding is the same as 
in Fig. [6] 



The LS 5039 lightcurve and spectra were modelled using a 
simple-minded leptonic model. The electrons are assumed to be 
accelerated efficiently in a small zone in the vicinity of the com- 
pact object with a standard y^ p power-law. Radiative losses due 
to inverse Compton emission and synchrotron emission gener- 
ate a distinctive steady-state electron distribution in this environ- 
ment dominated by stellar photons. The distribution has a promi- 



nent hardening between the energy at which inverse Compton 
losses enter the Klein-Nishina regime (yxN ~ 6 10 4 in LS 5039) 
and the energy at which synchrotron losses take over (ys ~ 10 7 
for a 1 G field). This is for instance the distribution found in the 
vicinity of the pulsar wind shock but it applies equally well to 
any leptonic model where particles are accelerated close to the 
compact object. The magnetic field was allowed to vary as the 
inverse of the orbital separation, as expected from a pulsar wind 
nebula. The model has only three parameters: the intensity of the 
magnetic field, the normalization of the electron distribution and 
the slope p of the injected power-law y e p . 

The cutoff in the very high energy gamma-ray spectrum is 
very sensitive to the magnetic field intensity, via the location of 
7s in the electron distribution. Fitting the high-state spectrum 
seen by HESS gives a rather constrained magnetic field inten- 
sity at periastron of 0.8+0.2 G. This value compares well with 
the values found using simple pulsar wind models which give 
5 {E-aO'iY^Rii G, where is the pulsar spindown power in 
units of 10 36 erg s 1 , crj, is the ratio of magnetic to kinetic energy 
in the pulsar wind in units of 10~ 3 and Rn is the distance of the 
shock to the pulsar in units of 10 11 cm. Fitting the HESS high- 
state spectrum also sets the injection slope to p = 2+0.3, close to 
the canonical value for shock acceleration. The normalization of 
the electron distribution implies an injection rate of 10 36 erg s _1 
for a radiative zone of 3 10 11 cm. These results are remarkably 
consistent with the expectations for a pulsar wind model. 

The spectrum is also found to fit extremely well the EGRET 
observations, adding credence to the reliability of this simple 
approach. The model predicts a strong variation in the GLAST 
band with a softening from high to low flux below a GeV (where 
synchrotron emission dominates the spectrum) but a hardening 
above a GeV (where inverse Compton emission dominates the 
spectrum). The HESS low-state spectrum is not explained to sat- 
isfaction. The model fits nicely the EGRET measurements but 
produces too many gamma-rays at 5-10 TeV. A possible solu- 
tion is a more complex orbital phase-dependence of the electron 
distribution at selected phases. Another solution is that the low- 
state spectrum corresponds to phases of strong attenuation and 
that emission from the created pairs contribute significantly to 
the spectrum. Additional HESS observations near minimum flux 
would be welcomed. 

The orbital modulation of the HESS emission is easily re- 
produced. A well-defined peak is predicted between phases 0.7- 
0.9 for which evidence may already be seen in the data. The 
lightcurve at GLAST energies is anti-correlated with the HESS 
lightcurve and has a peak at periastron, where the stellar photon 
density is maximum, and a minimum at inferior conjunction be- 
cause of the anisotropic effects in inverse Compton scattering. 
The GLAST spectrum below 1 GeV should be influenced by the 
tail of the synchrotron emission from the highest energy elec- 
trons. The peak synchrotron emission is at about 100 MeV for 
maximally accelerated electrons, regardless of magnetic field. 
Hence, if this component is detected, it will provide evidence 
that electrons are indeed accelerated with extreme efficiency in 
this source. 

Similar results for the magnetic field intensity and particle 
energy are found when a lower inclination is used, i.e. imply- 
ing a black hole compact object rather than a neutron star. In 
this case, the emission is thought to arise from a relativistic jet 
powered by accretion onto the black hole. Within the assump- 
tions of this work on the particle distribution, it is difficult to 
argue that a significant part of the emission occurs far along a jet 
since this does not naturally reproduce neither the spectrum nor 
the lightcurve measured by HESS. Most of the emission should 
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still occur close to the compact object. However, unlike in the 
case of a pulsar wind nebula, there is no independent theoreti- 
cal expectations in support of the magnetic field intensity (cer- 
tainly smaller than its equipartition value in the accretion flow) 
and particle energy that are derived. Therefore, the pulsar wind 
nebula model appears favoured independently of other possible 
considerations. 

Despite the complexity of the phenomena involved in pulsar 
wind nebula emission, it is found that the peculiar environment 
of a gamma-ray binary, most prominently the enormous lumi- 
nosity of the massive companion, severely constrains the num- 
ber of degrees-of-freedom in the model. A simple model suffices 
to reproduce most of the observations. The value of the magnetic 
field at the shock is found to be tightly constrained by the HESS 
observations to 0.8+0.2 G and the injection spectrum slope to 
p = 2 ± 0.3. These results confirm that gamma-ray binaries are 
promising sources to study the environment of pulsars on very 
small scales. 
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Appendix A: Inverse Compton spectrum for a mono-energetic beam of photons 

The purpose of this Appendix is first to carry out the integration set out in Eq. © and second to give an expression valid in the 
Klein-Nishina regime for the total spectrum emitted by a single electron scattering a mono-energetic beam of photons (Eq.|7]). The 
fraction of scattered photons per time, energy and steradian is given by Eq. ©, which can be expanded using Eqs. (Q][5]l 
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where primed (unprimed) quantities are measured in the electron (observer) frame, /jq> = cos ©' = /i'/ij + sin ff sin Q' x cos(0j - <p'), 
// = cos 6', yUo = cos 6q, fi' = cos ff Q etc. Re-arranging the last Dirac and performing the three integrations yields 
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The integration over Oi to obtain the full spectrum of radiation emitted by the electron is simplified if y e » 1. In that case, 
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(A.3) 



which is equivalent to saying the outgoing photon is emitted along the direction of electron motion when y e » 1 ■ The last Dirac can 
then be rewritten as a function of /ii : 
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with x = 



The integration over Q.\ is now straightforward, giving for the total spectrum: 
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Relativistic kinematics gives the domain of variation of the scattered photon energy e\ in the observer frame. The maximum e + and 
minimum e_ energies in the spectrum are : 
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The angle dependence of the maximum energy in the Thomson regime is (1 -Phq). For high electron energies, in the Klein-Nishina 
regime, the maximum photon energy is limited to y e m e c 2 and becomes almost independent of angle. 
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